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Abstract. The nebular properties (electronic density, ex- 
tinction, age, 0/H abundances) of Hii regions found along 
the bars of the sample of barred spiral galaxies studied by 
Martin & Friedli (1997) are examined. From line ratio di- 
agnostic diagrams, it is showed that regions located along 
the major axis of the bars have a normal photoionization 
spectrum, that is, line ratios reproductible from nebular 
conditions and ionizing star radiation field normally en- 
countered in extragalactic H ii regions. There is an indi- 
cation, however, that their degree of ionization might be 
somewhat different. Another ionization mechanism (high- 
velocity shocks or hard UV radiation) is clearly present 
for regions found nearby the centers of the galaxies. The 
electronic density of the regions along the bars is very 
close to that of disc regions (( N^) ~ 80 cm~^). On aver- 
age, bar and disc regions have a similar visual extinction 
{Av ^ Imag) with exceptions for some regions located 
near the bar dust lanes of the earlier types of galaxies in 
our sample. Although the average Ha equivalent width 
of bar Hll regions (^250 A) is half that of disc regions, 
this disparity could be due to uncertainties in the galactic 
continuum and does not translate into a significant age 
difference. The oxygen abundance distribution was also 
investigated in the bar of these galaxies. The 0/H scat- 
ter was found to be very small (<0.1dex) indicating that 
mixing of the chemical composition by gas flows is very 
efficient in a barred structure. 

Key words: Galaxies: abundances - Galaxies: ISM - 
Galaxies: spiral - Galaxies: starburst - Galaxies: kinemat- 
ics and dynamics - ISM: Hii regions 
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1. Introduction 

Galactic bars are the sites of highly diversified star forma- 
tion activities. Phillips (1993, 1996), Garcfa-Barreto et al. 
(1996) and Martin & FriedU (1997, hereafter MF97) have 
all showed that along certain bars, mostly found in late- 
type spirals, star formation (SF) can be quite intense. For 
instance, the bar of the SBcd galaxy NGC 4731 has a total 
star formation rate (SFR) of about 1.5 Mq yr"! (MF97). 
However, in other cases, SF in the bar can be very weak 
or completely absent. This is the case for most of the bars 
in early-type barred spirals (e.g. NGC 1300, NGC 1512, 
NGC 3351). The origin of these differences is not yet com- 
pletely understood. Numerical simulations (Friedli & Benz 
1995, MF97, Martinet & Friedli 1997) and observations 
(Martin & Roy 1995) suggest that the existence of mas- 
sive star formation in certain bars is a relatively brief event 
(~0.5-1.0Gyr) in the evolution of a barred system and 
that it mostly takes place during the formation of the bar 
itself. The amplitude and duration of the event is, how- 
ever, controlled by a complicated combination of differ- 
ent physical processes. For example, the time-dependent 
bar evolution, initial gas content, amplitude of gas flows, 
and mechanical energy injected in the interstellar medium 
(ISM) by supernovae ejecta are all factors that can influ- 
ence the level of SF activity in bars (MF97). Hence, it is 
essential to acquire more data on the properties of the H II 
regions formed in such environments to better constrain 
the relative importance of these factors, and consequently 
significantly improve SF recipes used in numerical simu- 
lations. 

The morphology, SFRs and other properties of SF 
along the bars of a sample of eleven spiral galaxies were 
studied in the first paper in this series (MF97). Includ- 
ing the large diversity in SF activity, MF97 found that 
the distribution of H ii regions can be highly asymmetri- 
cal in the bar. For some SBc spirals, large regions can be 
present outside the bar major axis. This morphology sug- 
gests that the SF process along the bar of late-type spirals 
is a chaotic process, not strictly confined along the major 
stellar/gas orbits defined by the barred potential. In such 
a case, one could expect similar properties for these H ii 
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regions when compared to the disc star forming regions. 
The situation could be different for bars in earlier types of 
galaxies for which the star forming regions are generally 
located next to dust lanes (MF97) . For one galaxy of their 
sample (NGC7479), MF97 also estimated the amount of 
gas flowing in the bar and falling into the galaxy center. 
They found that possibly as much as 75% of the gas in the 
bar is not transformed into stars. Numerical simulations 
suggest that this number is very dependent on the pres- 
ence of bar-induced shocks in the star forming ISM. It is 
then critical to determine whether there is some signature 
of these shocks in the H II regions along the bars. 

In this paper, spectrophotometric data are used to de- 
rive the physical properties of a sub-sample of H ii regions 
found along the bars of the galaxies studied by MF97. As 
a comparison sample, a few regions located in the discs 
of these galaxies have also been analysed. Using different 
diagnostic line ratio diagrams as defined by Baldwin et al. 
(1981) and Veilleux & Osterbrock (1987), the excitation 
of regions located in bar environments is compared to that 
of "normal" disc regions (Sect. 3.1). A similar analysis is 
performed for the electronic density using the appropri- 
ate sulfur line ratio (Sect. 3.2). The distributions of the 
visual extinction for both populations are also studied in 
Sect. 3.3. Using the Ha equivalent width indicator (e.g. 
Leitherer et al. 1999), we also infer the approximate age 
of Hii regions (Sect. 3.4). Since large-scale mixing of the 
chemical composition by bar-driven gas radial flows oc- 
curs in galaxy discs (Martin & Roy 1994; Friedli et al. 
1994), the 0/H distribution in bar environments is also 
investigated (Sect. 3.5). 

2. Database 

2.1. Observations 

The long slit observations were conducted during two 
runs in April and October 1994 at the equivalent 4.5 
meter Multi-Mirror Telescope on Mt Hopkins, Arizona. 
The "Blue Channel" spectrograph was used with a Lo- 
ral 3kxlk CCD. A grating with 500 grooves/mm and a 
blaze at 5410 A was employed; the spectral range covered 
in the reduced spectra is about 3500 A (-3500 ~7000) 
with a dispersion of 1.17 A/pixel. An order sorting filter 
was used (UV-36) and all the data were obtained with the 
CCD binned by a factor of two in the spatial direction 
(0.6"/pixel). All these observations were carried out with 
an unvignetted 2"xl50" slit positioned close to the paral- 
lactic angle to avoid any light loss. Three exposures of 15 
or 20 minutes were obtained for each slit position in the 
bar. Between exposures, the alignment of the MMT six 
mirrors was verified and corrected, if necessary. Numer- 
ous standard stars were observed with a larger slit (5") 
during the night for the flux calibration procedure. 

Table 1 presents the journal of observations for the 
H II regions of ten of the eleven objects studied by MF97 



(NGC5068 was not observed). For some objects, a few 
slit positions were required to optimize the number of H ii 
regions observed in the bar and the disc. Average seeing 
during these observations was about 1.2" and conditions 
were photometric. Note that the only galaxy in the sample 
considered as an AGN (LINER) is NGC7479. 

2.2. Data reduction and analysis 

The long slit spectra were reduced following standard pro- 
cedures available in the longslit package in iraf|^. First, 
a bias and flat-field correction was applied. The illumina- 
tion pattern along the slit was corrected using a set of sky 
fiat-fields taken in the same optical configuration. Wave- 
length calibration was done using a Helium-Neon- Argon 
exposure taken immediately after the science observation. 
Geometric distortion and alignment of the spectra were 
corrected by 2D-mapping of the spectral lines from the 
calibration sources using the fitcoords and transform 
algorithms. 

The next step, sky subtraction, was done by extract- 
ing a sky background using sections along the long slit 
outside the galaxy. To make sure that the signal from the 
galaxy was minimized, the average sky backgrounds were 
compared between the different slit positions since most 
of the slits used for the disc H ii regions were less contam- 
inated by the disc emission. The spectra were then flux 
calibrated using a sensitivity function obtained from a set 
of standard stars observed through the nights. Extraction 
of the spectrum for each Hii region was done using the 
spatial proflle at Ha seen along slit. In general, no contin- 
uum trace was found for these H li regions. The extraction 
trace was forced across the spectral range using the posi- 
tions of the main nebular lines along the spectral domain. 
Finally, the individual spectra were combined. 

The integrated fluxes of the main nebular lines H/3, 
[Oii]A3727, [Qui] AA4959, 5007, [Oi]A6300, [Nii]A6584, 
Ha, and [Sii] AA6717,6731 were measured using a Gaus- 
sian fitting algorithm available with SPLOT in iraf. A 
background continuum estimated from each side of the 
lines was automatically subtracted. 

Many H ii regions in bars show a strong underlying 
Balmer absorption at H/3. McCall et al. (1985) have shown 
that adding about 2 A of equivalent width for normal disc 
regions constitutes an appropriate correction. This correc- 
tion was applied for all the H/3 fluxes for the H ii regions 
in our control sample and the bar regions with shallow 
underlying absorption. However, for about 20% of the bar 
sample, this correction was not sufficient. For these re- 
gions, mostly located in the bars with a strong continuum 
(e.g. NGC 3504, NGC 5921, NGC 7479), the amplitude of 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation 
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Table 1. Journal of observations. 



GalaxA' 


Epoch 


# Slits 


Slit PA 


Exp()sur(> Tii!K>s s] 


<, Aimia.s.s> 


Ab„r" 


Adi-.-" 


L(,„,. [kpc] 


NGC 1073 


1994 Oct 9 


1 


68 


3 X 1200 


1.2 


2 


2 


5.0 


NGC 1087 


1994 Oct 9 


1 


139 


3 X 1200 


1.4 


2 


2 


1.6 


NGC 3319 


1994 Apr 17 


1 


40 


1 X 1200 


1.0 


2 


3 


4.0 


NGC 3359 


1994 Apr 16 


1 


26 


3 X 1200 


1.2 


8 


5 


5.8 


NGC 3504 


1994 Apr 16 


2 


140; 175 


2 X 3 X 900 


1.1 


6 


7 


6.8 


NGC 4731 


1994 Apr 16 


1 


126 


3 X 900 


1.3 


6 


4 


15.2 


NGC 4900 


1994 Apr 16 


1 


145 


3 X 900 


1.3 


4 


6 


3.4 


NGC 5921 


1994 Apr 16 


1 


159 


3 X 1200 


1.2 


5 


4 


8.1 


NGC 7479 


1994 Oct 4 


2 


2; 11 


2 X 3 X 1200 


1.3 


7 


9 


16.6 


NGC 7741 


1994 Oct 4 


2 


92; 98 


2 X 3 X 1200 


1.1 


7 


1 


3.5 



" Number of bar and disc Hll regions with detected H/3. 
Total length of the bar based on distance and measurements given in MF97. 



the underlying absorption was estimated from a Gaussian 
fit. In general, we found that about 5 A of equivalent width 
were necessary to assure a good correction. This value was 
applied for all bar regions with absorption higher than 
normal. 

All of the line fluxes were corrected for interstellar red- 
dening by comparing the Ha/H/? ratio to the theoretical 
Balmer decrement (2.86) for Case B recombination (that 
is, for nebulae with large optical depths for the Hi reso- 
nance lines) at 10* K. In reality, the temperature of bar 
Hii regions is probably aromid 7000 8000 K due to their 
high 0/H abundances (see below) so that our extinction 
based on the Balmer decrement might be overestimated by 
about 0.1 mag. Since the Balmer ratio represents only an 
approximation of the real extinction, we did not take this 
difference into account. The reddening law formulated by 
Savage & Mathis (1979) was assumed for reddening cor- 
rection. 

It is generally difficult to evaluate the accuracy of the 
absolute spc;ctrophotonietric fluxes for individual objects. 
The main uncertainties include the contamination by the 
bright galactic continuum, the underlying absorption at 
H/3, the accuracy of the photometric calibration, the spec- 
trum extraction in a crowded field of H ii regions, the posi- 
tioning of the slit, the correction for the interstellar extinc- 
tion and, of course, the intensity of the lines. We evaluate 
the accuracy of the fluxes to be at about 20-30%. Since we 
will only compare the global behavior of H II regions, these 
high uncertainties should not influence our conclusions. 

3. Nebular properties 

3.1. Excitation 

The basic source of ionization in an H ii region is the UV 
radiation field from young massive stars. The emission 
line spectra resulting from a pure photoionization field in 
a Stromgren sphere can be studied using diagnostic dia- 
grams (Baldwin et al. 1981; Evans & Dopita 1985; Veilleux 
& Osterbrock 1987; Osterbrock 1989; Dopita & Suther- 
land 1995; Rola et al. 1997). These standard diagrams 



are based on nebular line ratios like [O iii]/H/3, [O ij/Ha, 
[Nii]/Hq, and [Sii]/Ha. For high redshift galaxies, the 
red part of the spectrum is shifted to the near- infrared. 
If optical studies are needed, other diagnostic lines like 
[O ii] A3727 or [Ne ill] A3869 can also be relied upon (Rola 
et al. 1997). These diagnostic diagrams are extremely use- 
ful in distinguishing between normal photoionized regions 
and regions with another ionization mechanism (e.g. high- 
velocity shocks, hard UV fields). The line ratios above are 
also independent of the reddening correction and depend 
only on the accuracy achieved for the spectrophotometry. 
However, despite the usefulness of these diagnostic dia- 
grams for determining whether another ionization mech- 
anism is present or not, it is very difficult to identify this 
mechanism. For this, sophisticated nebular models must 
be used (e.g. Stasihska 1990; Dopita & Sutherland 1995). 

Figure 1 illustrates the first diagnostic diagram with a 

correlation between [Oin]/H/? and [Nil] /Ha for the bar 
regions (dark symbols) and the disc regions of our control 
sample (open symbols). The dotted lines define the limits 
of the sequence of normal Hii regions that can be foimd 
in Osterbrock (1989) and Kennicutt et al. (1989). The full 
line is the separation between normal photoionized regions 
and regions with another ionization mechanism. We have 
also indicated the effect of high-velocity shocks and mag- 
netic fields on these line ratios from models by Dopita & 
Sutherland (1995). The average location of Seyfert 2 and 
LINER galaxies in this diagram are also displayed. At this 
time, it is still not entirely clear that the spectral charac- 
teristics (mainly the high [Nil] /Ha ratio) of these last ob- 
jects are due to high UV radiation or high-velocity shocks 
(see discussion by Dopita & Sutherland 1995). However, 
Fig. 1 shows that excepting for four regions, all the Hll re- 
gions in our sample are located inside the sequence of nor- 
mal Hll regions. These four regions are located close to the 
nucleus in NGC 3504 (a starburst galaxy) and NGC 7479 
(LINER). Thus, from this diagram alone, bar Hii regions 
do not exhibit any sign of high- velocity shocks or hard-UV 
radiation. Apart from these central peculiarities, there is 
no dependence on radius. 
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Fig. 1. Diagnostic diagram of our sample of bar and disc 
Hii regions based on Osterbrock (1989). The sequence of 
normal regions is represented by the dotted lines. The full 
line indicates the separation between normal regions and 
regions with another ionization mechanism. The squared 
and circular symbols correspond to the average location of 
Seyfert 2 and LINER galaxies from Rola et al. (1997). The 
arrow lines indicate the effect of high- velocity shocks and 
magnetic fields on the line ratios from models by Dopita 
& Sutherland (1995). From the beginning to the end of 
the arrow, the velocity for the shock varies from 150 to 
500kms~^. The magnetic field parameter is zero for the 
bottom arrow and 4 fiG cm^^^ for the top arrow 

Although [Nil] /Ha is a good diagnostic ratio for high- 
velocity shocks or very hard ultraviolet radiation (Veilleux 
& Osterbrock 1987; Dopita & Sutherland 1995), other ra- 
tios like [Sii]/Ha or [Oi]/Ha are more sensitive to these 
ionization mechanisms. Figure 2 illustrates another diag- 
nostic diagram: the correlation between [Oiii]/H/3 and 
[Sll]/Ha. The dashed and full lines are as in Fig. 1. Ar- 
rows also show the effect of the high- velocity shocks and 
magnetic fields. Although most of the Hii regions fall 
within the normal region sequence, there is a small num- 
ber of bar regions outside the sequence. Figure 3 shows 
the sequence between [Oiii]/H/3 and [Oi]/Ha. The lat- 
ter ratio was detected in about 67% of the sample of bar 
regions but only in about 40% of the disc region sample. 

From these diagrams, it is clear that most bar Hii 
regions do not exhibit any systematic evidence of high- 
velocity shocks (>150kms^^) or very hard UV radiation. 
Only circumnuclear regions exhibit obvious signs of an- 
other ionization mechanism (fast shocks and/or hard-UV 
radiation) , as expected from the work of Kennicutt et al. 
(1989). Nevertheless, if these conditions do not appear to 
be present in bar regions, Figs. 2 and 3 suggest that shocks 
with lower velocity or an abnormal UV photoionization 
field cannot be excluded. This possibility is discussed in 
Sect. 4. 
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Fig. 2. Diagnostic diagram comparing bar and disc Hii 
regions. The dotted lines are the sequence of normal re- 
gions defined by Osterbrock (1989) and the full hue sep- 
arates the normal regions from regions which are not ex- 
clusively photoionized. The square and circular symbols 
correspond to the average locations of the Seyfert 2 and 
LINER galaxies from Rola et al. (1997). The arrows are 
as in Fig. 1 except that the magnetic field for the top one 
is 2/xGcm3/2 
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Fig. 3. Diagnostic diagram of bar and disc Hii regions. 
The dotted lines indicate the sequence of normal H ii re- 
gions. The full line is the separation between normal and 
abnormally excited regions (Osterbrock 1989). The square 
and circular symbols correspond to the average location 
of Seyfert 2 and LINER galaxies from Rola et al. (1997). 
The arrows are as in Fig. 1 



3.2. Electronic density 

The electronic density in Hii regions can be accessed 
with the line ratio of the [Sii] doublet at 6717-6731 A, 
i.e. 7 = -''SVi? (•Qg|;gj.]-,j.Q(,]j 1989). The density can be 



derived from nebular models published by Blair & Kir- 
shner (1985). Since the calibration depends on the ncbu- 
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Fig. 4. Electronic density distributions for bar (top) and 
disc (bottom) Hii regions. The density conversion de- 
rived from the [S ii] doublet ratio is from Blair & Kirshner 
(1985). Vertical dotted lines indicate the mean values 




Fig. 5. Distributions of the interstellar extinction Ay de- 
rived from the Ha/H/3 line ratio for bar (top) and disc 
(bottom) Hii regions. Vertical dotted lines indicate the 
mean values 



lar temperature (TVe oc T^'^), we will assume T,, = 10'* K. 
This value is probably too high for the Hii regions in 
our sample with abundances higher than the solar value 
(sec Sect. 3.5). However, for comparison purposes, this 
approximation is appropriate. The electronic density dis- 
tributions for the bar and disc H ii regions are illustrated 
in Fig. 4. In both samples, several regions show line ratios 
that are very close or larger than the low-density limit 
(7 ^ 1.4). The electronic densities derived in this regime 
are very uncertain because the doublet ratio becomes only 
weakly dependent on for values larger than about 1.3. 

Kennicutt et al. (1989) found that nuclei Hii regions 
tend to possess higher electronic densities on average than 
disc regions, with both classes showing a large range of 
densities. As seen in Fig. 4, the case of Hii regions lo- 
cated in bars is different. No significant difference is ob- 
served between the distributions of electronic densities of 
both populations. On average, 7 w 1.33 ±0.02 (bar regions) 
and 7 ?a 1.31 ±0.03 (disc regions). Bar Hii regions do not 
show any compactness with respect to disc regions. 

3.3. Extinction 

As discussed in the Sect. 2.2, the visual interstellar ex- 
tinction of individual H II regions can be derived from the 
Ha/H/3 line ratio. The values derived, however, are ap- 
proximate since in reality, the real extinction is probably 

not distributed uniformly but is patchy. Figure 5 presents 
the distribution of the visual interstellar extinction from 



both the bar and disc H 11 regions. There is a considerable 
extinction in both populations of Hii regions. The mean 
values are Ay « 1.3 and Ay w 1.0 visual magnitudes 
for the bar and disc regions, respectively. This difference 
is probably not significant since the extinction derived in 
bar regions, based on the assumption that the nebular 
temperature is 10'' K, might be overestimated by about 
0.1 to 0.2 magnitude. 

The extinction values for the bar regions show a large 
dispersion. Most of the Hii regions contributing to the 
highest values are located in the bars of NGC 7479 and 
NGC5921. As noticed in MF97, these regions are located 
close to the strong dust lanes seen in these bars. Nev- 
ertheless, the overall behavior shows that the interstellar 
extinction as derived from the Balmer decrement is simi- 
lar for bar and disc regions. However, using IRAS obser- 
vations, Phillips (1993) has shown that for circumnuclear 
regions the extinction derived from the Ha/H/? ratio can 
be underestimated by as much as 2 magnitudes. The "uni- 
form screen" model assumed for the extinction is probably 
over-simplistic. Any line ratios (e.g. R23) or other quanti- 
tative properties (e.g. integrated fluxes) severely affected 
by the interstellar extinction should be interpreted with 
caution for regions located in the inner parts of galaxies. 

3.4. Ha equivalent widths 

The equivalent widths (EW) of the Balmer emission-lines 
provide a measure between the number of ionizing and 
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continuum photons emitted in the H II region. As such, 
the EWs depend strongly on the stage of evolution of the 
ionizing stars, the initial mass function (IMF), and the 
metahicity (Dottori 1981; McCall et al. 1985; Copetti et 
al. 1986; Bresolin & Kennicutt 1997; Bresolin et al. 1999; 
Leitherer et al. 1999). As shown by Copetti et al. (1986) 
and more recently by Leitherer et al. (1999), EW(Hq!) 
and EW(H/3) can both be used as age indicators for Hii 
regions. In disc galaxies, the distribution of EW(Ha) ex- 
tends from about 100 A to 1500 A with a median value 
around 400 A. No obvious correlation with the Hubble 
type is found (Bresolin & Kennicutt 1997). Assuming an 
instantaneous burst of star formation and a solar metal- 
licity, these values correspond to an age range between 
IMyr to 7Myr (Leitherer et al. 1999). 

The accuracy of the Balmer line EWs is mostly de- 
termined by the uncertainty in the level of the nebular 
continuum which is severely contaminated by the galactic 
continuum. Because we could directly measure the con- 
tribution from the galactic continuum on our "off" band 
images used in MF97, we have only measured the EWs 
for the Ha line. The EW(Hq!) is also less affected by the 
interstellar extinction and the underlying absorption. The 
fraction of the galaxy light contributing to the nebular 
continuum was estimated from two photometric apertures: 
one covering the integrated light of the H ii regions, and 
the other located on nearest area devoid of any Ha emis- 
sion (determined from the Ha images) . A correction factor 
was then applied to the EW values measured directly from 
the spectra. These correction factors vary from 1.1 to 20 
depending on the location of the H ii region. 

Figure 6 illustrates the EW(Ha) distributions of the 
bar and disc Hll regions. The mean values for the distri- 
butions differ by about a factor of two: EW(Ha) « 250 A 
(bar regions) and EW(Hq;) w 560 A (disc regions). Fol- 
lowing the models of Leitherer et al. (1999) for an instan- 
taneous burst of star formation with a Miller-Scalo mass 
function and solar metallicity, the mean age of bar regions 
is about 5.3 Myr while disc regions are about 4.0 Myr old. 
However, the difference in age is less (< 1 x lO'^yr) when 
a Salpeter function is used to describe the IMF. Also, no 
age gradient seems to exist along the sequence of bar HII 
regions. This is an indication that HII regions should be 
ignified all along the bar and not only at bar ends with a 
subsequent migration towards the center. 

In their study of nuclear H ii regions, Kennicutt et al. 
(1989) found that the EW(Hq;) of the Hii region nuclei, 
with a median value around 25 A, are approximately 20 
times lower on average than that of normal disc regions. 
Such a difference cannot easily be explained by assuming 
that the correction for galactic continuum was underesti- 
mated. The authors rather favor the idea that the stellar 
continuum is high due to continuous star formation in the 
same region or an unusual stellar mass spectrum in the 
ionizing clusters. In the present case, however, it is diffi- 
cult to completely discard the effect of the contamination 
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Fig. 6. Distribution of the Ha equivalent widths EW(Ha) 
for the bar (top) and disc (bottom) H ii regions. The scale 
in the bottom panel indicates the approximate age of the 
Hii regions from models with a Salpeter IMF, instanta- 
neous burst of star formation and solar metallicity (see 
Leitherer et al. 1999 for details). Vertical dotted hues in- 
dicate the mean values 



of the galactic continuum to the nebular continuum to 
explain the discrepancy observed between bar and disc 
regions. The location of the aperture used to measure the 
galaxy continuum has a strong effect on the correction 
performed. The light distribution in bars is not uniform 
and some bars have very patchy dust features. A system- 
atic error of a factor 2 cannot be ruled out. In any case, 
no firm conclusion based on the difference observed in the 
EW(Ha) of both populations of H ii regions can be drawn 
from the actual sample. 

3.5. 0/H abundance (within bars) 

The oxygen abundance in Hii regions can be derived ei- 
ther through semi-empirical calibrations (e.g. Edmunds & 
Pagel 1984; McGaugh 1991; Pagel 1997) or directly when 
the temperature can be measured from the nebular lines 
[Oi]A4363 or [Nii]A5755. The latter, however, are gen- 
erally detectable only for H ii regions with low oxygen 
abundance. In our case, almost all the regions have solar 
or above-solar oxygen abundances; semi-empirical tech- 
niques have to be used. Even if the uncertainties related 
to these methods are generally quite large (±0.2 dex), it 
is worthwhile to derive the 0/H values to address the im- 
portant question of mixing of the ISM in bars. It is now 
well established that bars induce large-scale mixing of the 
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chemical composition in the disc of spirals (Martin & Roy 
1994; Zaritsky et al. 1994; Friedh et al. 1994). The radial 
flows of gas formed by bars flatten the strong (negative) 
abundance gradients generally observed in unbarred late- 
type disc galaxies. The importance of the homogenization 
effect is related to the bar strength as shown by Martin 
& Roy (1994). Very strong gas flows {v ^ 100 km s^^) are 
taking place along bars; efficient mixing should be also ob- 
served and the O/H scatter between the bar Hii regions 
should be smaller than what is observed in normal galaxy 
discs (See Sect. 4). 

The oxygen abundances for our sample of bar H ii re- 
gions were determined using three line ratios: [N ii]/[0 ill] , 

tr\ 1 /TTO J n [O II] A3727+[0 III] AA4959, 5007 mi 

[Oiii]/H/3, and R23 = ^ ^ 4ir^^ ■ 

conversion to relative oxygen abundances was done us- 
ing the calibration of Edmunds & Pagel (1984). Much has 
been written on the accuracy of these line ratios as abun- 
dance indicators (e.g. McGaugh 1991; Martin & Roy 1995; 
Stasinska 1998). For our sample, Fig. 7 compares the dif- 
ferent O/H values derived with all three indicators. The 
O/H values derived from [Nii]/[Oiii] are slightly higher 
(~0.1dex) than the values derived from [Olll]/H/3 and 
R23 for 12 -I- log(0/H) < 8.9. For 12 + log(0/H) > 8.9, 
that is, the abundance regime for most of the bars in our 
sample, the abundances from [Olli]/H/3 and [Nli]/[Oili] 
are slightly below the values given by i?23- These results 
were previously discussed by Martin & Roy (1994) and 



are due to discrepancies in the semi-empirical calibrations. 
For our purposes, we use the O/H values derived from the 
[Nii]/[Oiii] indicator; our conclusions are not affected by 
this choice. 

The distribution of the oxygen abundance along the 
bars of our sample is illustrated in Fig. 8. It is clear that 
the O/H scatter observed in these bars is well smaller 
than ±0.1 dex or even less. Martin & Belley (1996, 1997) 
have shown that the azimuthal O/H dispersion observed in 
the discs of normal and barred galaxies is generally equal 
or larger than ±0.2 dex. These abundance variations are 
probably not intrinsic but are in fact a combination of real 
inhomogeneities and uncertainties associated with empir- 
ical techniques (Roy & Kunth 1995). However, a compar- 
ative analysis remains possible when the same method is 
used to derive the O/H values (Martin & Belley 1997). 
Clearly, the chemical composition of H 11 regions in a bar 
is well-homogenized. This indicates that an efficient mix- 
ing of the chemical composition is taking place in the bar 
region (see next section). 

4. Discussion 

Degree of ionization. As discussed in Sect. 3.1, bar Hii 
region spectra do not exhibit any obvious signs of high- 
velocity shocks or hard UV radiation. However, there is 
marginal evidence from Figs. 2 and 3 that the ionization 
might be different for some bar regions. The degree of 
ionization at a specific position in a nebula can be accessed 
through the ionization parameter: 



U = 



Q(g°) 
Anr^cNp 



(1) 



where Q(H^) is the number of ionizing photons per 
unit time by the central source, r is the position in 
the nebula, and c the speed of light (Osterbrock 1989). 
The most obvious evidence that bar regions are differ- 
ent from disc regions is seen in the [O i] A6300 line. As 
described in Evans & Dopita (1985), the [O i] line is emit- 
ted in the transition zone of an Hii region which con- 
tains a significant fraction of neutral hydrogen. The line 
is then stronger when U is lower (or the ionizing stel- 
lar temperature is lower). Since we have detected the 
[O i] line in a much larger fraction of bar regions than 
in disc regions, this suggests that the ionization param- 
eter could be lower in the former population. Figure 9 
shows the correlation between [O i] A6300/[O ill] A5007 
and [On] A3727/[Olll] A5007 which is particularly depen- 
dent on U (Evans & Dopita 1985). A correction for the in- 
terstellar extinction has been applied to these line ratios. 
The bulk of bar Hil regions is located at ?7~0.0005. For 
the disc regions, the scatter is quite large but on average 
?7~0.001, larger than the value for bar regions. Unfortu- 
nately, our sample is not large enough to firmly confirm 
that U is indeed different for both populations of H 11 re- 
gions. 
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Fig. 8. 0/H abundance distributions along the bars of the galaxies in our sample. The 0/H values were derived using 

the [Nll]/[Olll] line ratio. The horizontal axis is the position of the regions normalized to the bar radius given in 
MF97. The open and black symbols differentiate the regions with regard to what side of the bar they are located. The 
abundance interval is 0.06 dex for galaxies displayed in the left column and 0.2 dex for those in the right column 



If a difference in the ionization parameter is really 
present, this could be due to many factors: differences 
in the initial mass function, age, richness of the OB as- 
sociations, or spatial distribution of the ionized material 
(Evans & Dopita 1985). As recently shown by Rozas et al. 
(1999), the luminosity function (LF) of the bar regions is 
much less regular than the LF of the disc regions in the 



strongly barred spiral NGC 7479. Their result, combined 
with our study on the nebular excitation, suggest strongly 
that the properties of the OB associations formed in bars 
differ from the normal associations of the disc. More work 
comparing LFs and nebular properties of a larger sample 
of bar and disc regions would allow us to investigate the 
origin of this difference. In the final paper in this series 
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Fig. 9. Diagnostic diagram showing the relation between 
the [Oi] A6300/[Oiii] A5007 and [On] A3727/[Oiii] A5007 
nebular line ratios. Both ratios have been corrected for 
the extinction. The curves represent models from Evans & 
Dopita (1985). The horizontal dashed lines show three val- 
ues for [/, 0.03, 0.003, and 0.0003. The vertical hues indi- 
cate the stellar temperatures used in the models, 50 000 K, 
40 000 K and 37 000 K. 

(Friedli & Martin, in preparation), we will also examine 
the properties of the clusters formed in diverse bar envi- 
ronments with high-spatial numerical simulations. 

Mixing and element production. The estimated timescale 
given by numerical simulations for which the star forma- 
tion activity phase lasts in bars is tsf 5 x 10^ yr (Mar- 
tin & Friedli 1997). How does this compare with mixing 
timescale? 

It is possible to roughly quantitatively evaluate the 
timescale of mixing of the ISM due to radial flows. Roy 
& Kunth (1995) have discussed the diverse mixing mecha- 
nisms of the oxygen abundance in the ISM in galaxy discs. 
Assuming a pure radial mixing due to gas flows funnelled 
in the bar, the upper limit for the time for gas to diffuse 
a length scale, Axrad, in the radial direction is: 



rad 



vl 



(2) 



where v is the radial flow velocity, and I is the mean free 
path for molecular clouds. In a typical bar, Aa;rad = 5kpc 
(see Table 1) and i; = 100kms~^. The value of the mean 
free path for the gas clouds is not a well-defined quantity 
in bars. In galaxy discs, I — 300 — 1000 pc (Roberts & 



Hausman 1984; Roy & Kunth 1995). If we assume / = 
500 pc for bars, we find Tj-ad ~ 5 x 10^ yr. However, since 
radial flows in bars are not stationary, the real mixing 
timescale could be even of the order of Trad = Axrad/'^, i-e. 
~ 5 X 10''yr. Putting all this together yields the following 
reasonable interval for the mixing timescale: 5 x 10^ ^ 
Trad 5 x 10® yr. Thus, Trad is shortcr than tsf meaning 
that the abundance content in H ii regions formed during 
this phase must be homogenized. 

It is also instructive to make rough (i.e. close-box) es- 
timates of the global abundance increase during Tgp as 
well as of the abundance fluctuations in H ii regions which 
should result from their age spread. If is the initial 
mean gaseous abundance in the bar region, e the global 
star formation efficiency, and Yz the net yield for the 
species considered, then the final mean abundance is given 
by: 



Zf = Z' + e(l - e)-^Yz 



(3) 



Interestingly enough, Z^ does not depend on the initial 
gas mass fraction. For instance, for the oxygen with an 
yield Yq ~ 0.006, an initial solar abundance « 0.01, 
and a typical SF efficiency e « 0.25, then Z^' w 0.012. The 
global increase of oxygen abundance is thus only about 
0.1 dex. Equation 3 can in fact also be applied to each in- 
dividual H II regions with exactly the same numbers; the 
fluctuations in the oxygen abundance are thus expected 
to be of the order 0.1 dex, which is indeed what is ob- 
served (Sect. 3.5). However, we do not observe any clear 
trend between the age and metallicity for bar H ii regions. 
In dwarf galaxies, the metal enrichment of Hii is not ob- 
served and metals are probably locked in the hot phase, 
(see e.g. Tenorio-Tagle 1996; Kobulnicky 1998). The situ- 
ation could be similar for star forming regions in bars. 

5. Summary 

The main results concerning the properties of H ii regions 
located within bars can be expressed as follows: 

1 ) From standard diagnostic diagrams, the excitation of 
most Hll regions appears normal and similar to the one 
of disc regions. There are some exceptions for nuclear re- 
gions where an ionization mechanism other than photoion- 
ization seems to be present. However, there is marginal 
evidence that the ionization parameter in bar regions is 
lower than in disc regions, suggesting that the properties 
of the OB associations might be different in bar environ- 
ments. 

2) The electronic density distribution as derived from the 
[S ii] line ratio is similar to that observed for normal disc 
regions. The mean density is ~ 80 cm^'^. Star forma- 
tion regions in bars have the same "compactness" as disc 
regions. 

3) The Ha/Hp extinction indicator reveals that, on aver- 
age, bar regions have a visual extinction Ay ^ l.^mag, 
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0.3 mag more than disc regions. This difference is mainly 
due to the fact that some regions are located near the bar 
dust lanes of the earlier types of galaxies in our sample. 

4) The average Ha equivalent width for the bar regions 
is about 250 A, half that of disc regions. While this could 
indicate an older population for bar Hll regions, the cor- 
responding age difference is probably too small to be sig- 
nificant owing to the large uncertainties introduced by the 
galactic continuum correction. 

5) The 0/H abundance distribution of these Hu regions is 
remarkably hom,ogeneous. This is the result of the gaseous 
radial flows in bars inducing mixing of the ISM, as seen 
on a larger scale in the discs of barred spirals. 
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